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Abstract 

<D 

Motivated by the Exact Parity Model and other theories, the hypothesis that each of 
the known neutrinos oscillates maximally with a sterile partner has been put forward as an 
explanation of the atmospheric and solar neutrino anomalies. We provide detailed predic- 
tions for muon and electron flux ratios induced in the Kamiokande and SuperKamiokande 
detectors by sub-GeV atmospheric neutrinos. Several different, carefully chosen cuts on 
momentum and zenith angle are proposed, emphasizing the role of up-down flux asymme- 
tries. 



The solar |]T] and atmospheric [|] neutrino observations and the LSND || experiment 
provide strong evidence that neutrinos have nonzero masses and oscillate. Many specific 
neutrino mass scenarios have been proposed as explanations for some or all of these results. 
In this paper we will focus on the hypothesis that each of the known neutrinos oscillates 
maximally with an exactly or effectively sterile partner. Though non-minimal, this scenario 
is theoretically well motivated by the Exact Parity Model |J and other theories ||. 

Denote the three sterile neutrinos by v' e , u' , and v' T . In the Exact Parity Model, all known 
particles are paired via an exact parity symmetry with mirror partners in order to achieve 
a parity symmetric dynamics. Since parity eigenstates must also be energy eigenstates, 
the mass eigenstate neutrinos are maximal mixtures of u a and v' a (where a — e, //, r) in 
the limit of small intergenerational mixing. The models of Ref. provide approximate 
maximal mixing via an approximate discrete symmetry of the neutrino mass sector. We 
will assume small intergenerational mixing, analogous to the quark sector, for the purposes 
of this paper. 

Maximal v e — v' e and — v' oscillations have been put forward as explanations of 
the solar and atmospheric neutrino anomalies respectively [|, ||. This scenario is also 
compatible with the LSND experiment. Maximal v e — v' e oscillations can solve the solar 
neutrino problem for 8m 2 ee , in the large range [|7|, §, 

3 x 1(T 10 ~ \5m 2 ee ,\/eV 2 ~ 0.9 x 1(T 3 , (1) 

where the upper bound is the most recent experimental limit . Maximal — v' oscilla- 



tions can explain the atmospheric neutrino anomaly |T0[ provided that [El 11 



lO- 3 ^\6m^,\/eV 2 ^10-\ (2) 

Note however that the atmospheric neutrino experiments are sensitive to both — v' 
and v e — v' oscillations in principle. The implications of this for atmospheric neutrino 



experiments were discussed qualitatively in Ref . |12| , where rather crude analytic estimates 
were provided of some relevant measurable quantities. The purpose of this paper is to 
provide a detailed quantitative study of the implications of maximal — v' and v e — v' e 
oscillations for the Kamiokande and SuperKamiokande experiments. 

We will restrict our analysis to sub-GeV neutrinos. The main reason for this is as 
follows. Consider the zenith-angle averaged "ratio of ratios" (R), where 

obs / o \ 

(N,/N e ) MC {6) 

The quantities N ejfl are the numbers of electron- and muon-like events. The numerator de- 
notes numbers actually observed, while the denominator the numbers expected on the basis 
of a Monte-Carlo simulation without oscillations. Since (R) is measured to be significantly 
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less than 1, we conclude that the v e oscillation length must be much larger than the v ^ os- 
cillation length. Given this, v e oscillations are expected to impact more significantly on the 
sub-GeV data than the multi-GeV data, because the neutrino oscillation length decreases 
with energy. We will not comment further on the multi-GeV neutrinos, except to note that 
the Kamiokande and especially the preliminary SuperKamiokande data fl3| are consistent 



with maximal — u' oscillations and are inconsistent with the minimal standard model. 

In the water- Cerenkov Kamiokande and SuperKamiokande experiments, neutrinos are 
detected via the charged leptons £ a (£ a — e or //) produced primarily from quasi-elastic 
neutrino scattering off nucleons in the water molecules: v a N — > £ a N' (a = e, /i). The total 
number N(£ a ) of charged leptons of type £ a produced through u a N — > £ a N' is given by 

r^C 1 P Q ni ax r i r +i 



POO r Q ni ax 

N(£ a ) — Ut / dE / dq d cos i/j / dcos9 
d 2 F a (E J 9) d 2 a a {E,q, cosip) 



dE dcosO dq dcosip 



P(u a ^u a] E,9). (4) 



Here d 2 F a /dEd cos 9 is the differential flux of atmospheric neutrinos of type v a of energy 
E at zenith angle 9. The term is the effective number of target nucleons. The function 
d 2 a a /dqd cos ip is the differential cross section for quasi-elastic scattering, v a N — ► £ a N', 
where ip is the scattering angle relative to the velocity vector of the incident v a (the az- 
imuthal angle having been integrated over), and q is the energy of the charged lepton £ a . 
The function P[y a — > v a ; E, 9) is the survival probability for a v a with energy E after trav- 



eling a distance L = \J (R+ h) 2 — R 2 sin 2 9 — R cos 9, where R is the radius of the Earth and 
h ~ 15 km is the altitude at which atmospheric neutrinos are produced. It can be obtained 
by solving the Schrodinger equation for neutrino evolution including matter effects (for a 
review see, for example, Ref . [ ) . In our analysis we have a pair of two-flavour oscillation 



subsystems, v e — v' e and v^ — v'^ where the vacuum mixing angles are maximal (sin 2 2#q 6 ' = 1 

and sin 2 2^q M =1). There are thus only two free parameters: Sm 2 e , and Sm^,. Note that 
Eq.(f|) must be modified in order to obtain zenith angle and momentum binned charged 
lepton events. 

For the case of sub-GeV neutrinos, almost all the information necessary for the right 
hand side of Eq.@ is available in the published references. We have used the differential 
cross section d 2 a a /dqd cos ip of Ref.|T5|. The differential flux of atmospheric neutrinos 



d 2 F a / dEd cos 9 without geomagnetic effects is given in |L6|], but we have used the differential 
flux which includes geomagnetic effects JT7|]. 



For reference and for completeness, we first briefly discuss the case where only maximal 
— v' oscillations occur. Figure 1 shows the variation of (R) for sub-GeV neutrinos with 
<5m 2 ^,. The preliminary SuperKamiokande central value plus the la confidence level band 

is also plotted. The datum implies 10~ 3 ~ bvn 2 ^ j eV 2 ~ 10 -2 at the la level, with a much 
larger range allowed at the 2a level (but also constrained by the zenith angle dependence of 



2 



the multi-GeV sample). Note that while a major thrust of this paper is to examine possible 
manifestations of v e — v' e oscillations, it is quite possible that 5m 2 e , is too small to lead 
to any observable effects for atmospheric neutrinos. In that case, — v' oscillations will 
be the only effect occurring. This case is almost indistinguishable from the oft-considered 
maximal z/ M — v T oscillation solution to the atmospheric neutrino anomaly [fill. Note that 



future long baseline z/ T -appearance experiments should help distinguish the two cases. 



We now turn on maximal v e — v' e oscillations. As discussed qualitatively in Ref.[|12j, there 
are several interesting effects that can occur when both v e — v' e and — u' oscillations are 
present. One simple effect of v e — v' e oscillations, for a given 5m 2 r , is to increase (R) 
relative to what it would be in the absence of these oscillations. We illustrate this in 
Fig. 2, which gives (R) as a function of 5m 2 e , for various values of Sm 2 ^,. We have used 
the usual Kamiokande momentum cuts 0.1 < p e /GeV < 1.33 and 0.2 < p^/GeV < 1.5. 
The current experimental measurements of (R) by the Kamiokande and SuperKamiokande 
collaborations are 11131 , 



(R)Kam = 0.6018$ ± 0-05 

(R)sKam = 0.635 ± 0.033 ± 0.053 (preliminary). (5) 

In the absence of v e — u' e oscillations (or equivalently, when \5m 2 e i | ~ 2xl0~ 5 eV 2 ) the above 
measurements suggest |5m^„/| ~ 5 x 10~ 3 eV 2 . However, as Fig. 2 shows, much larger values 
of Sm 2 i can fit these data equally well when v e — v' e oscillations also occur. For example, an 
(R) value of about 0.6 is obtained for 5m 2 , ~ 6 x 10~ 2 eV 2 when 5m 2 e , ~ 10 -4 eV 2 . Figure 
2 also shows that all of the range of Eq.(|l|) is consistent with the (Super)Kamiokande data. 

The observation that v e — v' e oscillations increase (R) is interesting for future long 
baseline neutrino oscillation experiments. If u e — v' e oscillations are numerically important 
for atmospheric neutrinos, then the long baseline determination of Sm^, should not be 
consistent with that implied by atmospheric neutrinos under the assumption that — v' 
oscillations only are occurring. A discrepancy would signal that an additional effect is at 
work for atmospheric neutrinos, with v e — v' e oscillations being a candidate. 

One way to experimentally uncover the presence of both — v' and v e — v' e oscillations 
is through the energy dependence of (R) . It turns out that a more sensitive way is through 
the zenith-angle dependence of fluxes |12[ , as we now explain. 

Consider the quantities Y" e ° and Y® where 

y o = (N-/N+)\ osc o _ (N-/N+)U 

e " (iV-/iV+)|„ " (N-/N+)\ no _ osc - 1 ] 

Here N~ (AT+) is the number of electrons produced in the detector with zenith angle 
cosO < (cosO > 0). iV^ are the analogous quantities for muons. (The zenith angle 
O for the charged leptons should not be confused with the zenith angle 9 defined earlier 
for neutrinos.) The numerators are the predictions for N~/N+ (a = e,fi) in the model 
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while the denominators are the same quantities in the absence of oscillations. (Note that 
N~ I 'N+\ no _ osc are close to 1 for symmetry reasons, but not exactly equal to it.) The 
numerators and denominators, being ratios of fluxes, should be approximately free of the 
systematic errors arising from uncertain absolute fluxes (±30%), uncertain \x to e flux ratios 
(± a few percent) and the uncertain cross-section. The Y® are a measure of up-down flux 
asymmetries. It is important to observe that Y® depends only on Sm 2 e ,, and Y^ depends 
only on 5m 2 m ,. 

We have plotted F e ° as a function of 5m 2 e , in Fig. 3 We have used the usual 



Kamiokande cut 0.1 < p e /GeV < 1.33 as well an the alternative 0.5 < p e /GeV < 1.33 cut. 
In the former case we see that Y" e ° decreases from 1 to about 0.82 as 5m 2 e , is increased from 
10 -5 eV 2 towards the experimental limit. The effect is much more pronounced when the 
alternative cut on momentum is made, with Y" e ° decreasing to about 0.67. The reason for 
this is the greater correlation between the zenith angle distribution of the electrons and the 
zenith angle distribution of the incident u e flux when the average energy of the sample is 
increased. A price has to be paid in statistics since the neutrino flux is larger at smaller 
energies, but a gain is made in the size of the signal relative to the underlying systematic 
error (which is at the few percent level for flux ratios). This is important because the 
statistical error will be reduced as the data sample grows, whereas the systematic error 
will remain the same (barring great improvements in the theoretical prediction of neutrino 
fluxes) . 

Figures 4 and 5 are similar to Fig. 3, except we have redefined iV+ to be the number 
of electrons with cos0 greater than 0.2 and 0.6, respectively. Similarly N~ is redefined to 
be the number of electrons with cosG less than —0.2 and —0.6, respectively. We use the 
notation F e 2 and Y® s for the quantities corresponding to Eq.|| Note that these cuts on the 
zenith angle also increase the magnitude of the effect, again with a concomitant decrease 
in statistics but a gain in signal relative to systematic error. 

Figures 6-8 are the corresponding figures for the muon case, Y^. Note that in this case 
the usual Kamiokande momentum cut is 0.2 < p^/GeV < 1.5 and we have considered the 
alternative cut of 0.5 < p^/GeV < 1.5 also. The peak at high Sm 2 ^, is due to an oscillation 
node in the downward-going neutrinos. For the muon case there is already interesting 
evidence for a zenith angle dependence of the flux of multi-GeV neutrinos. In this context, 
we would like to emphasize the important role that the ratios Y^ and Y e , with the various 
momentum and zenith angle cuts, could play in determining whether the atmospheric 
neutrino problem is due to oscillations alone or to both z/ M and v e oscillations. Current 
data clearly show that oscillations are the dominant effect. However, the presence of non- 
trivial v e oscillations cannot yet be ruled out and will be further tested as SuperKamiokande 
continues to acquire data. 

Note the insensitivity of Y^ to Sm^, for the quite large range 3 x 10~ 4 ~ Sm^/eV 2 ~ 
10~ 2 . The corresponding values of about 0.75 — 0.80 and 0.65 — 0.70 for Y® are, in a sense, 
the most "likely" ones. This insensitivity may provide a useful test of maximal v„ — v' 
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oscillations and may also be useful in discriminating against other possible explanations of 
the atmospheric neutrino anomaly. 

We will now quantify the decrease in statistical error expected with time. After about 
1 year of running, the preliminary SuperKamiokande data show statistical errors for Y e ^ 
of about 8% and 12% for the standard momentum cut and alternative momentum cut 
lower limit of 0.5 GeV, respectively. These should decrease to about 3 — 4% and 5 — 7%, 
respectively, after 5 years of running. These are small enough to see or rule out many of 
the effects plotted in Figs. 3-8. 

In conclusion, we have provided a quantitative analysis of the implications of the hy- 
pothesis that both the and v e neutrinos are maximally mixed with sterile partners for the 
(Super)Kamiokande atmospheric neutrino experiments. Maximal mixing, motivated by the 
Exact Parity Model and other theories |J , implies that the oscillations can be described 
by just two parameters, 5m 2 e , and 8m 2 ,. As Figs. 2-5 show, if |5j7ig e /| ~2x 1CT 5 eV 2 , then 
the v e — v' e oscillation length is too long to affect the experiments. As is well known, the 
experimental results can then be explained by maximal — v' oscillations with 5m 2 > in 
a range around 5 x 10~ 3 eV 2 . Zenith angle asymmetries for the detected sub-GeV muons 
produced are then expected (see Figs. 6-8). (Significant zenith angle dependence for multi- 
GeV muons has, of course, already been observed.) For \8m 2 e , | ~ 2xl0" 5 el/ 2 , the effects of 
u e — u' e oscillations become significant. They have two main effects. First, they can increase 
(R) (see Fig. 2) relative to what it would be in their absence. Second, they can lead to zenith 
angle asymmetries for the electrons produced in the detector (see Figs. 3-5). While the cur- 
rent data are consistent with the minimal case of maximal z/ M — v' oscillations only, the 

possibility that the electron neutrino also oscillates maximally with \5m 2 e , I ~ 2 x 10" 5 eV 2 
is an interesting possibility. Furthermore, this possibility will be tested stringently in the 
near future as more data is collected, especially at SuperKamiokande. 

Note added 

While this paper was in the final stages of preparation, a paper by J. W. Flanagan, J. 
G. Learned and S. Pakvasa appeared on the hep-ph archive (9709438), which also discusses 
the importance of up-down zenith angle asymmetries of the charged leptons induced by 
atmospheric neutrinos. 
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Figure Captions 

Figure 1. The zenith angled averaged "ratio of ratios", denoted (R), plotted as a function 
of Sm 2 ^, for the case where v e — v' e oscillations are unimportant (5m 2 e , ~ 2 x 1CT 5 eV 2 ). 
The usual SuperKamiokande momentum cuts have been assumed. The band denotes the 
lcr allowed region from preliminary SuperKamiokande data. 

Figure 2. (R) plotted as a function of 5m 2 e , for several values of 5m 2 , with the usual 
momentum cut used by SuperKamiokande to define its sub-GeV sample. The solid lines, 
going from top to bottom, correspond to Sm^/eV 2 = 1(T 3 , 3 x 1(T 3 , 6 x 1(T 3 , 1(T 2 , 6 x 1(T 2 
and 3 x 10~ 2 , respectively. The band denotes the la allowed region from preliminary 
SuperKamiokande data. 

Figure 3. The quantity F e °, defined in the text, as a function of 5m 2 e ,. The solid line uses 
the usual cut 0.1 < p e /GeV < 1.33 on the electron momentum p e , while the dashed line 
uses the alternative cut 0.5 < p e /GeV < 1.33. 

Figure 4. Similar to Fig. 3 except we have plotted F e 2 defined in the text. 

Figure 5. Similar to Fig. 3 except we have plotted F e a6 defined in the text. 

Figure 6. The quantity Y®, defined in the text, as a function of bm 2 ^,. The solid line uses 
the usual cut 0.2 < p^/GeV < 1.5 on the muon momentum p^, while the dashed line uses 
the alternative cut 0.5 < p^/GeV < 1.5. 

Figure 7. Similar to Fig. 6 except we have plotted Y^ 2 defined in the text. 
Figure 8. Similar to Fig. 6 except we have plotted F° 6 defined in the text. 



